QPOs and firehose instabilities in neutron star magnetospheres in accreting systems 
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We show that the interaction of an accretion disk with the magnetosphere of a neutron star 
can excite resonant shear Alfven waves with Hz-kHz frequencies in a region of enhanced density 
gradients. This is the the region where accretion material flows along the magnetic field lines in 
the magnetosphere. We argue that due to the pressure anisotropy produced by the plasma flow, 
firehose instabilities are likely to occur. Furthermore, for a dipolar field topology, we show that a 
new instability develops due to both magnetic field curvature and the plasma flow. 
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Recent observations of quasi-periodic oscillations 
(QPOs), particularly kHz QPOs, in X-ray emissions from 
accreting binaries, have aroused a lot of interest in the 
astrophysical community. The QPOs are very strong and 
remarkably coherent with frequencies ranging from ~ 10 
Hz to ~ 1200 Hz. They have been observed in the X-ray 
flux of about 20 accreting neutron star sources and five 
black hole sources by Rossi X-Ray Timing Explorer. Al- 
most all sources have shown twin spectral peaks in the 
QPOs in the kHz part of the X-ray spectrum, with the 
value of the peak separation being anti-correlated with 
the QPO frequencies [1-3]. 

The clear similarities of kHz QPO properties in black 
hole systems to those in neutron star binaries [4, 5], and 
in white dwarf systems [6], suggest that the oscillation 
likely originate in the accretion flow surrounding the cen- 
tral object. Motivated by this argument, a variety of 
accretion-based models have been proposed. Most mod- 
els, however, fail to provide a general explanation of QPO 
features in all potential sources. See, for example, Li & 
Narayan [7] and references therein. 

Li & Narayan [7] studied Rayleigh- Taylor and Kelvin- 
Helmholtz instabilities at a possible interface between the 
star's magnetosphere and the accretion disk. They found 
that modes with low order azimuthal wavenumbers are 
expected to grow to large amplitude and to contribute 
to kHz QPOs. In their study, the magnetic field al- 
lows an interface with abrupt spatial discontinuities in 
the flow density and/or angular velocity across the in- 
terface. They ignore, however, the dynamical role of the 
magnetic field of the central object and its interaction 
with the accretion disk/flow. 

Based on theoretical models and observations of the 
aurora in the Earth's magnetosphere, Rezania et al. [8] 
and Rezania & Samson [9] (hereafter paper I) have re- 
cently proposed a generic magnetospheric model for ac- 
cretion disk-neutron star systems to address the occur- 



rence and the behavior of the observed QPOs in those 
systems. In the Earth's magnetosphere, the occurrence 
of aurora is a result of the resonant coupling between 
shear Alfven waves and fast compressional waves (pro- 
duced by the solar wind). These resonances are known 
as field line resonances (FLRs). Paper I argued that this 
resonant coupling is also likely to occur in neutron star 
magnetospheres, due to interaction with the accreting 
plasma. The MHD interaction of the infalling plasma 
(with a sonic/supersonic speed) with the neutron star 
magnetosphere, would alter not only the plasma flow 
toward the surface of the star, as assumed by current 
QPO models, but also the structure of the star's mag- 
netosphere. The magnetic field of the neutron star is 
distorted inward by the infalling plasma of the Keplerian 
accretion flow. Furthermore, the plasma would likely be 
able to penetrate through magnetic field lines and pro- 
duce enhanced density regions within the magnetosphere 
(similar to the magnetospheric interface introduced in 
[7]). Any instability in the compressional action of the 
accretion flow would alter the quasi-equilibrium pressure 
balance between the inward pressure of the infalling flow 



pvf/2 and the outward magnetic pressure 



Bli 



*Also at Institute for Advanced Studies in Basic Sciences, Zan- 
jan 45195, Iran.; Electronic address: vrezaniaOphys .ualberta. ca; 
URL: http: //fermi .phys .ualberta. ca/~vrezania 



This process then can excite perturbations in the en- 
hanced density region. Here p and v r are the density 
and radial velocity of the infalling matter and B p is the 
poloidal magnetic field on the plane of the disk. In anal- 
ogy with the interaction of the solar wind with the earth's 
magnetosphere [10], one would expect the excitation of 
resonant shear Alfven waves, or FLRs, [9]. 

Paper I assumed a simple geometry with a rectilin- 
ear magnetic field, and showed that, in the presence of 
a plasma flow, two resonant MHD modes with frequen- 
cies in the kHz range can occur within a few stellar radii. 
The results in paper I gave a reasonable prediction, both 
quantitative and qualitative, of the kHz oscillations ob- 
served in the X-ray fluxes in X-ray binaries. The exis- 
tence of a non-zero plasma displacement along the mag- 
netic field lines, which oscillates with resonance frequen- 
cies and then modulates the flow of the plasma toward 
the surface of the neutron star, can explain the kHz quasi- 
periodicity in the observed X-ray flux from the star. See 
[9] for more details. 
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In this letter we examine the above problem in a more 
realistic configuration for the star's magnctosphere. We 
base are model on the fact that the topology of the mag- 
netospheres of isolated pulsars is likely close to dipolar. 
In accreting pulsars, however, the geometry of the mag- 
netosphere will be distorted due to the inward flowing 
plasma, specifically on the plane of the accretion disk. 
Nevertheless, we approximate the topology of the mag- 
netosphere of an accreting neutron stars with a dipolar 
geometry in order to study the structure of shear Alfven 
waves in the presence of an ambient flow. 

Furthermore, due to the existence of ambient flow 
along the magnetic field, the plasma pressure is not ex- 
pected to be isotropic. This can be understood by noting 
that the plasma pressure parallel to the field lines can be 
defined as 



P\\ 



m 



J f(u\\ - v p ) 2 d 3 u, 



(1) 



where U| = u • B/B is the thermal velocity of the plasma 
particles and v p is the ambient flow velocity along the 
field lines. Here /(r, u, t) is the plasma distribution func- 
tion satisfying the Vlasov equation [11]. Similarly, the 
perpendicular component of the pressure can be calcu- 
lated from 



Pi- 



rn 



j fu 2 ± d 3 u/2. 



(2) 



Hence, the scalar pressure of isotropic MHD must be re- 
placed by a diagonal pressure tensor with two compo- 
nents: a parallel component p\ | acting along the field lines 
and a perpendicular component p± acting in the perpen- 
dicular direction. The latter can be considered to be the 
ram pressure. Consequently, in the MHD equations, the 
flow pressure must be written as: p = p±l+ (p± —p\ )bb, 
where I is the identity tensor and b = B/B is the unit 
vector along the magnetic field line. 

The linearized, perturbed magnctohydrodynamic 
equations with anisotropic pressure in the presence of an 
ambient flow now have parallel and perpendicular com- 
ponents: 



/cMv 



V dt 



dSv 
~dt 
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v • V(5v + 5v ■ Vv ) =-V||5p|| 

II 

V||5B N 



+(P± -P||) 
v • V<5v + <5v • Vv 



B 

V_l<5 tp ± 



„ (SB VB + B SVB 



V x (5v x B + v x SB), 



dt 

V • SB = 0, 

5p\\ = -2p\\SB\\/B, 
Sp±_ =p±SBu/B, 



(3a) 
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(3b) 

(3c) 

(3d) 
(3e) 
(3f) 



where <5v = d£ /dt, 5 = 1 + 2 (c j_ — c 2 | ) jv\ , and we ignore 
the perturbation in the plasma density, i.e. Sp ~ 0. Here 
P>P|hPJ-! v i an d B are the unperturbed quantities while 
Sp\\, Sp±, <5v, and <5B are the perturbed quantities. Equa- 
tions (3e) and (3f) are calculated from the two equations 
of state for pi | and p± that are known as double adiabatic 
equations, f t (p l{ B 2 / p 3 ) = and £ t (p ± /(pB)) - [12]. 
S = 1 if the pressure is isotropic, i.e p\\ = p± = p. 



To avoid complexities, we shall ignore the rotation of 
the star and consequently neglect both the toroidal field 
B^ and velocity v$. These assumptions simplify our cal- 
culations significantly, and, we believe, allow the model 
to retain the important physics. Furthermore, we do not 
consider a jump condition in the enhanced density re- 
gion as discussed by [7]. As a result, we do not address 
Rayleigh- Taylor and/or Kelvin-Hcmholtz instabilities. 



We expand the MHD equations (3) in the orthogo- 
nal coordinate system (p,v,(j)), where p — cosd/r 2 is 
the magnetic field-aligned coordinate, v = sin 2 9/r nu- 
merates magnetic shells in the direction perpendicular 
to the field line, and <f> the is azimuthal coordinate. 
The components of the metric are = h^h^, h u — 
r 2 j sin Q\J 1 + 3 cos 2 9, and h<$, = rsind where (r,9,(f>) is 
the spherical coordinate. The metric component al- 
lows a convenient representation of the dipolar magnetic 
field in the form B = /i mag //i M e M where // mag is the mag- 
netic dipole moment of the star. Note that, for a dipolar 
magnetic field B = B p e p , B p h p = /i mag is constant which 
leads to V x B = 0. We further assume that the ambient 
flow velocity is along magnetic field lines, ie. v = i> p e M . 



Assuming 



5{p, v, 4>, t) <~ 5(v)e 



ikfj, —tut imtfi 



(4) 



we can reduce the equations of motion (3) to one second 
order differential equation for Sv u as 
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d 2 Sv L , . .d5v v 



G{p,v)5v v = 0, 



F{p,v) = 



1 



G<ji,u) 



(4 + v 2 A )rji 
1 

(4 + 



(ci + v\)d„{Tn +m) + 2E-^r] 1 d I/ ln^ + 



2ikv p r]i(3c^ — c\)d v hih h 
h^iuD + - V • v p ) 



/^Su^ fc 2 + (^lnfe, y ) 2 
ft, 2 «Wfl + AT„ - V • v p 



+ — - =— ( tA(3C|| - c ± )i/2 - d^h^/K J 

»7i = {htQ/hJ/ftiwD - KJQ - (m 2 /h 2 )(cl + v\) {iw D + K 4> -V- v p )], 



■j-dvihfjitt,) - ^-d u h^ + -i— {h^d u v p - Vpd u h^) - 



h v 



h^h v 
2 ~ 1 



iu) D + K v - V • v p 



Q = lo 2 d + K 2 + (icon - K^V ■ v p - -j^[k 2 + -^{d^f}, 



(5) 
(5a) 



(5b) 
(5c) 
(5d) 

(5e) 



where wd = w- kv p /h^ is the Dopplcr shifted frequency, 
(i = (i, v, <f>), va = B/y/Anp is the Alfvcn 



Ki = v p • V In h, 
wave velocity, C| = y/p\\/p and c± = y/p±/p are sound 
velocities parallel and perpendicular to the direction of 
the magnetic field. We note that in deriving Eqs. (5) we 
assumed that d v \p± + B 2 / (8ir)] ~ 0. 

The shear Alfvcn resonance happens at rji = or 
equivalently at Q = leading to 

v 2 " 1 

Lu 2 D +iV-Vp u D +Kl-K 4> V-Vp—^f[k 2 +-^{d ti h 4> ) 2 ] = 0. 



(6) 



As a result, resonance frequencies will be given by 

w± = top//i M - (i/2)V • vp ± (1/2)VK, (7) 
A = 4t£H[fc 2 + (0„ ln/i ) 2 ]//i 2 - (V • vp - 2tf*) 2 . 

For a rectilinear configuration, resonance cigcnfrcqucn- 
cies Eq. (7) for an incompressible plasma flow along 
the field lines, V • v p = 0, with an isotropic pressure, 
P\\ = p± = p, reduce to ones we obtained in paper I, and 
as expected: u>± — k(v p ±va)- 

Equation (7), however, shows that whenever the am- 
bient flow is compressible, i.e. V • v p ^ 0, and/or A < 0, 
waves do not propagate and an instability develops. In 
general, the ambient plasma is fairly incompressible, i.e. 
V • Vp = 0. However, due to the topological deforma- 
tion of the magnetosphere caused by the compressional 
action of accreting material, a non-zero density gradient 
through the plasma, and so V ■ v p ^ can be expected: 



1 dp 

V ' v p = — 17 

p at 



-v 



p 



ldp 
~ ~~pdi 
Via p. 
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7 P 



V Pl 



(8) 



As a result, the resonant mode will grow (decay) if v p • 
Vp > (< 0). Approximating the plasma inflow velocity 



with the free fall velocity v p ~ v s (r) = (2GM/r) 1 / 2 , the 
growth/decay timescale will be as order of r ~ r/v p = 
(r 3 /2GM)V2 ^ 6 x 1CT 5 (M/M )- 1/2 (r/lOkm)^ 2 s. 
Therefore, the closer to the star the faster the instability 
develops. 

The condition A < is satisfied whether 



I : 

II : 



H < 
(V-\ 



2 ^ 2 
Cm > C 



vl/2, 



2K^f >Av 2 A \E\[k 2 



(9a) 

(0 M ln^) 2 ]/j£. 

(9b) 



Case I, that is known as the firehose instability in the 
literature, happens when p\\ is much larger than p± +pu, 
where pm — B 2 /(8ir) is the magnetic pressure. 

The magnetic field channeling the parallel plasma 
streams experiences a similar instability. Whenever the 
flux tube is slightly bent, the flowing plasma exerts a 
centrifugal force, that tends to enhance the initial bend- 
ing. The field line bending is proportional to the den- 
sity of energy in plasma motion along the magnetic field 
<~ pv 2 ~ Recalling Eq. (1), the ambient plasma flow 
would enhance the firehose instability in the magneto- 
sphere of an accreting neutron star: 



II 



' 2 , 2 
C||+V 



(10) 



where c' 2 | 
Maxwcllia] 

ishes. Inserting Eq. (10) into Eq. (9a), we find 



— (raj ' p) J fu 2 d 3 u and we assume that a 
Maxwellian distribution function, so the cross term van- 



v 2 > v 2 j2 + {cl-c'l). 



(11) 



Therefore, a firehose instability develops whenever con- 
dition (11) is satisfied. 

For an isotropic pressure, i.e. S ~ 1, the firehose in- 
stability would not be expected. However, an instability 
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may arise if the condition (9b) is satisfied. For an incom- 
pressible flow, we find that Eq. (9b) reduces to 

v p > Wl + ? 2 , (12) 

where q = k / (d ^ In h^). It is necessary to note that 
this instability is enhanced whenever both v p ^ and 

ln/i^ 7^ 0. The latter condition is due to the curvature 
of magnetic field lines. Therefore, the non-flat topology 
of the magnetic field can trigger some MHD instabilities 
through the magnetosphere. An interesting note is that 
the condition (9b) is only valid for for an isotropic pres- 
sure flow. When the wave transfers energy to the flow, 
the wave decays, and the extraction of flow energy by the 
wave will lead to a growing mode. For a superAlfvenic 
flow, the extraction of energy from the flow and grow- 
ing MHD waves is very likely [13]. In this case also the 
growth/decay timescale will be as order of r ~ r/v p = 
(r 3 /2GMy/ 2 - 6 x 10" 5 {M /Mq)- 1 ' 2 ^ /lQkmf' 2 s. 

Furthermore, by approximating the Alfven velocity by 
va ~ B{r)/ 'yj '47T ■pi where ps = M/(vg 47rr 2 ) is the free 
fall mass density (v p ~ vs(r) = y/2GM/r), the instabil- 
ity condition (12) is satisfied for distance further than 

r > 1.7 x 10 6 cm (1 + q 2 ) 2 ' 7 $J M^ /7 {M/M Q )- 1 / 7 . 

(13) 

Therefore, this instability is very likely to develop at a 
position where the accretion disk can distort the dipolar 
magnetosphere. Here /i 2 6 is the magnetic field dipolc 



moment at the surface of star in units of 10 26 G cm 3 and 
Mi 7 is the mass of accretion rate in units of 10 17 g s _1 . 
We calculate Eq. (13) for a solar mass neutron star with 
10 km radius. 

To summarize, we believe we have shown that there 
are generic plasma instabilities associated with the ra- 
dial position where the inflowing material of the accre- 
tion disk leads to a distortion of the inner magnetosphere 
of the neutron star and field aligned plasma flows. We 
suspect that these instabilities may be relevant in the un- 
derstanding of some details of observed QPOs in X-ray 
binary systems, particularly when linked to shear Alfven 
waves. The MHD interaction of the infalling plasma with 
the neutron star's magnetosphere can alter the topology 
of the star's inner magnetoshere. The plasma flows and 
topology lead to the possibility of : (1) firehose insta- 
bilities associated the pressure anisotropy produced by 
the plasma flow; (2) convective growth of waves in the 
plasma flow. The unstable modes might produce shear 
Alfven resonances with large amplitude giving strong 
quasi-periodic variations in X-ray fluxes. 
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